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Abstract 

We describe the WA89 Collaboration experimental data on A, S + , H - , and Q~ 
baryons, and A and H antibaryons production in S~ collisions with C and Cu targets 
at 345 GeV/c (^s^n ~ 25.5 GeV) in the frame of the Quark-Gluon String Model. How 
the theoretical results compare to the experimental data is discussed. Finally, some 
relations among the values of the model parameters obtained with the help of quark 
combinatorics are presented. 



PACS. 25.75.Dw Particle and resonance production 
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1 Introduction 



The Quark-Gluon String Model (QGSM) and the Dual Parton Model (DPM) are based 
on the Dual Topological Unitarization (DTU) and they quantitatively describe [1-6] 
many features of high energy production processes, including the inclusive spectra of 
different secondary hadrons, their multiplicities and multiplicity distributions, etc., both 
in hadron-nucleon and hadron-nucleus collisions at fixed target energies. QGSM and 
DPM account for the main features of secondary production at collider energies. 

In the QGSM, high energy interactions are considered as proceeding via the exchange 
of one or several Pomerons, and all elastic and inelastic processes result from cutting 
through or between Pomerons [7J. Each cut Pomeron leads to the production of two 
strings of secondaries. The inclusive spectra of hadrons in the final state of the collision 
are related to the corresponding fragmentation functions of the quarks and diquarks 
at the end of the strings. These fragmentation functions are constructed by using the 
Reggeon counting rules [9]. 

To study the interaction with a nuclear target the Multiple Scattering Theory (Gribov- 
Glauber Theory) is used, and, thus, this interaction is considered as the superposition 
of the interactions of the incident hadron with different nucleons in the target [5j [8] . 

In previous papers [1-6], where the secondary production by proton, pion, and kaon 
beams was reasonably described, the model parameters were already fixed by comparison 
of the theoretical calculations with experimental data. 

Here we consider the production of secondaries in S~C and TrCu collisions and we 
present the comparison of the theoretical results with the experimental data at 345 GeV/c 
obtained by the WA89 Collaboration [10] . To include in our description interactions out 
of E~ beams we have to deduce the expressions both of the momentum distribution 
function of ds diquark in S~ hyperon, and of its fragmentation functions into secondary 
hadrons. 

However, in the case of S~ beam the agreement with the experimental data of the cal- 
culations obtained by using these standard fragmentation functions is not good enough. 
The agreement becomes better when some additional polynomial factors are included 
into our fragmentation functions, but it is not clear whether this is a result of some 
special structure of the strange baryons and their resonances, or it is simply connected 
to possible experimental inconsistencies. 

2 Production of secondaries on nuclear targets in 
QGSM 

In QCD hadrons are composite bound state configurations built up from the quark 
ipi(x), i = 1, ...N c , and gluon G%(x), a = 1, N% — 1, fields. In string models baryons are 
considered as configurations consisting of three strings attached to three valence quarks 
and connected at one point (small volume) x, called the "string junction" (SJ) [11-13]. 
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The corresponding wave function can be written as 



\B) = i/ji(xi) ■ if)j(x 2 ) ■ ^fc(x 3 ) • J l3k (x,xi,x 2 ,x 3 ) , 
J ijk (x,Xx,x 2 ,x 3 ) = %(x u x) ■ ^ 3 jl (x 2 ,x) ■ $%,(x 3 ,x) ■ e i,j ' k ' . (1) 

Here the operator x) represents the gluon field string with endpoints at X\ and x. 

Such a "star" (or Y configuration of the baryon wave function \B) is favoured [TTJ [13] 
with respect to the also possible "triangle" (or A) configuration. 

Let us discuss in more detail the processes in which one or several Pomerons are 
exchanged. Each exchanged Pomeron in hadron-nucleon and hadron-nucleus interaction 
corresponds to a cylindrical diagram that, when cut, produces two showers of secon- 
daries. The inclusive spectrum of secondaries is determined by the convolution of di- 
quark, valence quark, and sea quark distributions, u(x, n), in the incident particles with 
the fragmentation functions, G(z), of quarks and diquarks into the different hadrons. 
The diquark and quark distribution functions depend on the number of cut Pomerons, 
n, in the considered diagram. 

In the QGSM one calculates the invariant cross section 



xe da f E da 2 



a ine i dx J a ine i d 3 p 



d 2 PT , (2) 



where x = 2p\\/y/s is the Feynman variable xp, and xe = 2E/y/s, and one has then to 
use one value of (p%) (here we have taken the value (jjj) = 0.35 (GeV/c) 2 ) to make the 
transition to the values of da/dxp which are presented in the experimental papers [10] . 

Thus, for the case of a nucleon target the inclusive spectrum of a secondary hadron 
h in QGSM has the form [I]: 

' 1 ■ — = V w n ■ (j> h (x, n)+w D - <f> h D {x) . (3) 



Cjnei dx 71=1 

The functions <j) h (x, n) determine the contribution of diagrams with n cut Pomerons, and 
w n is the probability for this process with n cut Pomerons to occur |15j . The second 
term in the right-hand side of Eq. (3) describes the contribution of diffraction dissociation 
processes, where the triple- Reggeon diagrams are also included [H [3] . The expressions 
of (j) h (x, n) for TjN (N = p, n) collisions in Eq. (3) have the form [TJ IT3] : 

^ N (x, n) = n) ■ /*(x_, n) + f h q {x + , n) ■ /* (x_, n) + 2(n - l)^(x + , n) ■ / s ft (x_, n) , 

(4) 

where 



x± = -[y4my/s + x 2 ± x] , (5) 



with rriT = yw 2 + Pt the transverse mass of the produced hadron, and f qq , f q , and f s cor- 
responding to the contributions of diquarks, valence quarks, and sea quarks, respectively. 
The quantities fi(x + ,n) and fi(x-,n) account for the contributions to <p(x,n) of the S- 
hyperon beam and of the target nucleon, and they are determined by the convolution 
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of the diquark and quark distribution functions with the corresponding fragmentation 
functions, e.g., 

fg q (x+,n) = / u qq (xi,n) ■G^ q {x + /xx) dx x . (6) 

For the case of nuclear targets one has to consider the different possibilities of one 
or several Pomeron cuts in each of the v hadron-nucleon inelastic interaction blobs, as 
well as of cuts between Pomerons. For a YT A collision, one of the cut Pomerons links a 
valence diquark and a valence quark of the hyperon projectile with a valence quark and 
a diquark of one target nucleon, respectively, while the additional Pomerons link the sea 
quark- ant iquark pairs of the projectile, either with diquarks and valence quarks, or with 
sea quark-antiquark pairs, of the target nucleons. 

As one example, the diagram for the inelastic interaction of the S~-beam with two 
target nucleons is shown in Fig. [TJ In the blob of the £ _ iVi inelastic interaction one 
Pomeron is cut, while in the blob of the S~iV 2 interaction two Pomerons are cut. To 
include all diagrams, i.e. to account for all possible Pomeron configurations and per- 
mutations, is essential for a correct calculation. The process shown in Fig. [1] satisfies 
the condition that the absorptive part of the hadron-nucleus amplitude is determined by 
combinations of the absorptive parts of hadron-nucleon interactions, according to rules 
given in refs. [16-19]. 
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Figure 1: One of the diagrams for the inelastic interaction of an incident S -hyperon 
with two target nucleons N± and N2 in a S^A collision. 

For /iA collisions, where n inelastic interactions occur with v target nucleons, one has 
that n > v, and 1 < rti < n — v + 1, rij being the number of cut Pomerons connecting 
with the i-t/i target nucleon. By denoting the relative weight of the contribution with 
cut Pomerons in every hN blob as w^ N , and by using the same procedure as in ref. [5], 
we can write the corresponding expressions for the inclusive spectrum of the secondary 
hadron h produced in a collision, where all possible Pomeron permutations and all 
possible different quark contents of the protons and neutrons in the target have to be 
accounted for. 

In particular, the contribution to the inclusive spectrum of the diagram in Fig. [T] is 
written as follows: 
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xe da 



= 2Vg • wf* ■ wf N * ■ {f h qq {x + , 3) • (*_, 1) 

+ 3) ■ /£(*_, 1) + 3) ■ [/£(*_, 2) + 2) 

+ 2/f(x_,2)]} , 



(7) 



where is the probability of pure inelastic (non diffractive) interactions with v target 
nucleons of a nucleus A to occur. 

The diquark and quark distributions, as well as the fragmentation functions, are 
determined from Regge intercepts, and their expressions were presented in Appendix 1 
of ref. [20] (see also [2TJ [22] )• Now, for the case of the presence of a ds diquark in the 
beam, they are given in the Appendix of this paper (see below). 

For secondary baryon production, the diquark fragmentation function contains two 
contributions. The first one corresponds to the production from the sea of a BB pair in 
the midrapidity region (see Fig. |2]), and it will be discussed in detail in the Appendix. 



Figure 2: Cut-chain diagram corresponding to the diquark fragmentation function for 
the production of a central BB pair. 

The second contribution is connected with the direct fragmentation of the incident 
diquark into a secondary baryon with conservation of the SJ. In the frame of QGSM three 
possibilities exist for this second contribution [20J . The secondary baryon can consist of: 
(a) the SJ together with two valence and one sea quarks, (b) the SJ together with one 
valence and two sea quarks, and (c) the SJ together with three sea quarks. These three 
possibilities are shown in Fig. [3J 

The fraction of the energy of the incident baryon carried by the secondary baryon 
decreases from (a) to (c), whereas the mean rapidity gap between the incident and 
secondary baryon increases. 

The processes shown in Figs. and [3b are the standard ones in QGSM and DPM, 
and they determine the main contribution to the multiplicity of secondary baryons in 
the fragmentation region. 
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Figure 3: QGSM diagrams describing secondary baryon B production by diquark d: (a) 
initial SJ together with two valence quarks and one sea quark, (b) SJ together with one 
valence quark and two sea quarks, and (c) SJ together with three sea quarks. 

On the other hand, the diagram shown in Fig. [3b leads to the difference in baryon and 
antibaryon production at rapidities far from the incident baryon (baryon charge diffusion 
in rapidity space). The role of such a process in the description of experimental data 
was considered in detail in refs. [21-26]. For £v4 collisions at 345 GeV/c the relative 
contribution of this diagram is rather small. 

The fragmentation function of diquark d, with quark content d = qiq2, into a sec- 
ondary baryon B through the processes shown in Fig. [3k, Eb, and [3b has the form [20J: 



Here (3 = 1 — acsj, with as,j being the intercept of SJ Regge trajectory, e is the relative 
suppression factor of the (c) contribution with respect to the processes (a) and (b), and 
cle is a normalization parameter. In the present calculations we use the values e = 0.024, 
as j = 0.9, and = 1.33, as in ref. [21]. The factor (1 — z) 1 accounts for the fact that 
the intercept of the 0-meson Regge trajectory, pa 0, is smaller than the standard 
non-vacuum Reggeon intercept, ~ 0.5. The value of the parameter 7 is half the 
difference between the strangenesses of the considered diquark and that of the secondary 
baryon. The powers of z and (1 — z) are changed when either q\ or q2 is a strange 
quark. The values of vf for different quarks i and baryons B are determined by quark 
combinatorics [26| [27] . These values are presented in the Appendix (see below) . 

3 QGSM description of the experimental data 

The experimental data for A, £~, S + , H~, and fi~ baryons, A and S + antibaryons 
production in S~ collisions with C and Cu targets are presented in ref. [10J in terms of 
dcr/dxp. 

Let us start our analysis from the left panel of Fig. [U where the data on p spectra 
at pp collisions and their comparison with our old QGSM calculations are shown. The 
agreement is good, as it was usually obtained in our previous papers. 

For comparison, we also present similar predictions for the case Y,~p — > S~X, which 
is theoretically similar to pp — > pX, but about 2 times smaller in the region of moderate 



Gf J {z) = a N -z fi -K qa 'Z 



2.5-/3 




■(l-z)+v£.e.(l-z) 2 }-(l-zy . (8) 
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xp. It's generally accepted that this difference can be connected to the rather large 
probability to produce a A in the S~ case. At small xp one also has suppression of the 
£~ production by sea quarks, what leads to a larger difference in the two considered 
reactions. 




■ pp -> pX 400 GeV/c [29] 
* pp -> pX 1 75GeV/c [28] 
Tpp — ^ pX 1 0OGeV/c [28] 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 



10 r ■ r Cu -» r X 345 GeV/c [1 0] 
• 5T C r X 345 GeV/c[10l 
□ pCu-»pX p T =0.3GeV/c 100GeV/c[30] 
opC->pX p,=0.3 GeV/c 1 00 GeV/c [30] 
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Figure 4: Experimental Feynman-x distributions of secondary p produced in pp collisions 
at 100, 175, and 400 GeV/c, together with their description by QGSM, and the QGSM 
predictions for Y,~p — > H~X (dashed curve), left panel. Experimental Feynman-x distri- 
butions of secondary S~ produced in S~C and H~Cu interactions at 345 GeV/c [10], and 
similar distributions of secondary p on the same targets [30] , and their comparison with 
the the predictions for £~'s (solid curves) and for protons (dashed curves) by QGSM, 
right panel. 



The experimental data on secondary protons production on nuclear targets at 100 
GeV/c and pp = 0.3GeV/c [30] are also in good agreement with the QGSM, as it is 
shown in right panel of Fig. 4 and in [5]. On the contrary, experimental data for £~ 
production on C and Cu targets shown in the right panel of Fig. 0] are in contradiction 
with our calculations when we use the fragmentation functions directly taken from the 
Reggeon counting rules, as it was done in all previous papers [USEE]. One can 
immediately see that these data have different shape that the proton data for the same 
targets. Unfortunately, the proton data [30] presented here were measured at fixed 
Pp = 0.3GeV/c, and there are no proton production data on nuclear targets integrated 
over pp, which could be used for direct comparison. However, the difference in the shapes 
of the distributions seems to be too large, especially when keeping in mind that the data 
of [30] were successfully described by the QGSM in ref. [5]. 

Now, let us compare the spectra of secondary S~ and A obtained in ref. [ID] . The 
calculated ratio of these spectra, together with the corresponding experimental points 
obtained by using the data of ref. [ID] , are presented in Fig. 5. Both £~ and A should 
be produced in the interval Xp = 0.3-0.8 mainly by the process shown in Fig. [3^, when 
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the ds diquark fragments into them by picking up, either aim or a c/-quark from the sea, 
respectively. The fragmentation functions for both channels should be somehow similar, 
with differences coming from other contributions including triple-Reggeon terms. This 
comes from the fact that ds and dd diquarks appear in the incident S~ with probabilities 
2/3 and 1/3, respectively, and, on top of that, the fragmentation of the dd diquark into 
strange secondary baryons is suppressed by the strangeness suppression factor S/L (see 
Appendix). So, the ratio £~/A should be approximately constant in the considered 
interval. However, the experimental ratio of S~/A yields increases more than 3 times 
in the interval xf = 0.3-0.6, in total disagreement with our theoretical expectation. We 
want to stress that we have never met such a large disagreement with the experimental 
data in any of our previous calculations, and one has also to note that this disagreement 
with the experimental data on the ratio S~/A can not be corrected by any theoretically 
meaningful modification of the ds-diquark fragmentation functions. 




Figure 5: Experimental ratios of secondary £~ to A yields vs Xp in YrC and YrCu 
interactions at 345 GeV/c [10], together with the QGSM predictions for these ratio 
(solid curve). The dashed and dotted curves show the ratios when taking into account 
different contributions to the QGSM calculations (see details in the main text). 

The absolute value of S~ to A ratio in the QGSM is more model dependent. In 
the simplest approximation, when all final dds states are assumed to be S~ and all uds 
states are assumed to be A (secondary S° are usually registered as A after radiative 
S° — > A7 decay), the ratio S~/A is equal to 1. However, this can be changed by the 
resonance production, e.g. the state £*(1385), that can be produced with rather large 
cross section, has a dominant decay mode An that would transfer some part of dds states 
into additional A's, and it would thus decrease the S~/A ratio. 

To account for this effect we assume that when only ds diquark fragmentation is 
considered, the hyperon production leads to the following empirical rule for the £~/A 
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ratio [27]: fragmentation is considered, the S /A ratio would be. 

R ( * S J /A = XT/A = 0.6 . (9) 

Such an assumption has provided a reasonable description of the A xf spectra in pp, pre, 
kp collisions [201 ETJ [22] . Both E~ and A spectra should be affected by nuclear effects in 
a similar way. 

To show the structure of our calculations in more detail, in Fig. 5 we present the 
ratio of X - to A when only part of all contributions is accounted for. The dashed curve 
shows the result of calculations when only the first term in the right-hand side of Eq. (4) 
is considered. In this case the ratio of S~ to A in the fragmentation region (xf > 0.2) is 
very close to 0.6, the very small difference from this value coming from the contribution 
of dd diquark fragmentation. At small xp, the values of x_ in Eq. (5) increase, leading 
to the decrease of the functions /™(x_,n) in Eq. (4), and to the decrease of the £~/A 
ratio. 

The ratio of S~ to A production when only fragmentation of ds diquark is considered, 
and without diffraction dissociation nor triple- Reggeon contributions, is shown in Fig. 5 
by a dashed curve. If all diquark and quark terms in Eq. (4) are included (dotted 
curve), the ratio is smaller and the agreement with the experimental data becomes worse, 
indicating that the ds diquark is responsible for the value 0.6 of the S~ to A production 
ratio. The solid curve shows the result of the complete QGSM calculation with diffraction 
dissociation and triple- Reggeon contributions, and though now the theoretical curve goes 
up at large Xf, it still presents a very significant disagreement with the experimental data. 
It seems clear that the agreement of the QGSM predictions with the experimental data 
can not be obtained by simply considering a slight variation of any of the contributions 
in Eq. (4). 

As the QGSM can not reproduce the experimental ratio of X~ to A production, it 
can not describe separately the x^-distributions of both £~ and A, that present clearly 
different experimental behaviours. It is nevertheless interesting to find out which one, 
S~ or A distributions, or both, are at the origine of the disagreement in the ratio. The 
two xp distributions for S~ and A are presented in Fig. 6. 

The difference between the S~ and the A xf distributions could be explained by as- 
suming that in the case of A production, but not in that of S~ production, the resonance 
decay contribution is quite significant. In all cases, when including the resonance decay 
contribution for S _ production the xp distribution would become softer, i.e. narrower, 
and consequently the agreement with the experimental data would be worse. 

The quark and diquark distribution and fragmentation functions are given by the 
Reggeon counting rules. Thus, for some fragmentation function G qq (z) in Eq. (6) having 
asymptotical behaviours G qq (z — > 0) ~ z a and G qq (z — > I) ~ (1 — z) 13 , the fragmentation 
function is written as the simplest interpolation of these two asymptotical behaviours in 
the form [S]: 

G h qq (z) = a h -z a -(l- zf , (10) 

where is a parameter which determines the inclusive density of a produced hadron h. 
We will call such a form of the fragmentation function as the standard one. However, a 
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Figure 6: Experimental xp distributions of secondary £~'s (left panel) and A's (right 
panel), produced in S~C and TrCu interactions at 345 GeV/c [10J, together with the 
corresponding QGSM predictions. The solid curves show the result of the calculations 
with only the standard polynomial terms. The dashed curves in the left panel show the 
calculations with one additional polynomial factor in the third term of the ds fragmenta- 
tion function, and the dashed curves in the right panel show the result of the calculations 
when accounting for the resonance contribution to the spectra of A. 

slightly more complicate form with additional polynomial factors and new parameters b 
and b n is also possible: 

G h qq (z) = a h -z a -(l-zf-(l + b- z K ) . (11) 

The theoretical xp distributions of E~ and A shown by solid curves in Fig. 6 have 
been calculated by using the standard form of the QGSM fragmentation funactions, 
without any additional polynomial factor in Eq. (11). The distributions so obtained are 
in an evident disagreement with the data for the case S~ production, and the attempt 
to include an additional polynomial factor in the fragmentation function of ds diquark 
(dashed curves in the left panel of Fig. 6) does not lead to any significant improvement 
of the agreement with the data, since though the absolute values of the spectra increase 
the shapes remain being wrong. 

In the case of A production (right panels of Fig. 6), the agreement with the ex- 
perimental data of the QGSM calculation with the standard form of the fragmentation 
functions with the experimental data is reasonable, and it becomes better when taking 
into account that some part of the A-hyperons are produced after resonance decay and 
they consequently have smaller xp. To include this effect in our calculation we have in- 
troduced into the fragmentation function Eq. (1A.31I) ) the additional factor (1 — |). The 
QGSM results obtained with this modified fragmenation functions are shown in Fig. 6 
by dashed curves. 
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The spectra of secondary S + and H are shown on Fig. 7. These secondaries are 
produced with cross sections several times smaller than secondary E~ and A. 




Figure 7: Experimental xp distributions of secondary S + (left panel) and S~ (right 
panel), produced in YTC and in YTCu interactions at 345 GeV/c [ID] , together with the 
corresponding QGSM predictions. The solid curves show the result of the calculations 
with only the standard QGSM expressions. The dashed curves are obtained with with 
one additional polynomial factor in the fragmentation functions. 

In the standard approximation of the QGSM, for S + production only one valence 
s quark, Eq. (A. 14), as well as an s quark through the diagram in Fig 3b, from the 
incident £~ can be used to fragment into S + , but another possibility that is usually 
considered is the resonance production of E*°(1385) or A(1405) in the process of ds- 
diquark fragmentation (see Fig. 3a), and the subsequent decay into £ + 7r~. The results of 
the calculations in this approximation (shown by solid curves in the left panel of Fig. 7) 
are in reasonable agreement with the experimental data. The shapes of the curves 
are correct, though the normalization are underestimated on the level of 30%. This 
disagreement can be corrected (dashed curves in the left panel of Fig. 7) by introducing 
a polinomial factor (1 + 5-2) into the last term of the ds diquark fragmentation function 
in Eq. flA~35l) . 

For the case of S _ production (right panel of Fig. 7) the contribution of the ds- 
diquark fragmentation is decreased by the strangeness production suppression factor. 
The calculation of S~ production with only the standard terms in the diquark fragmen- 
tation functions (solid curves in the right panel of Fig. 7) results in a too fast decrease 
of the spectra when increasing xf- The calculation with an additional polinomial factor 
(1 + 3\fz) in the last term of the diquark fragmentation function Eq. (A. 38) leads to a 
better agreement with the experimental data (dashed curves on the right panel of Fig. 7), 
except for the region of low xf, where the model results are significantly higher than the 
experimental data. 
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In the case of secondary Q~ production by £~ beam, the incident s quark should 
pick up two strange quarks from the sea. The cross section of this process should clearly 
be small due to the presence of the squared strangeness suppression factor. The exper- 
imental points for secondary Q~ production are presented in Fig. 8. Here the standard 
QGSM predictions are in reasonable agreement with the data. 




Figure 8: The xp distributions of secondary fi~'s produced in S~C and Y^Cu interac- 
tions at 345 GeV/c [10J. The curves show the QGSM predictions with only the standard 
QGSM expressions in the diquark fragmentation functions. 

The experimental yields of A's and S +, s [TO], which only contain sea antiquarks and 
do not depend on the S J contribution, are shown in Fig. |9j The corresponding description 
by the standard QGSM fragmentation functions clearly underestimates the A yields at 
xp > 0.4. For the case of S + production the agreement is reasonable. The description 
of the A and H spectra obtained by using the fragmentation functions of eqs. (A. 20) 
and (A. 21), where an additional polinomial factor (1 + 20 • z 2 ) has been included, leads 
to a better agreement with the experimental data of the A spectra, letting apart some 
overestimation in the small xp region. 

4 Conclusions 

We present the QGSM description of the experimental data [10] on secondary hyperon 
production in E~C and TrCu collisions at 345 GeV/c. These data are of special interest 
because the main contribution to the spectra of secondary A, S - , and 5~ at xp > 0.3 
comes from the direct fragmentation Fig. 3a of the incident ds diquark, with rather small 
background from another subprocesses. 

From the invariant cross section in Eq. (3) and by using the value (pj) = 0.35 
(GeV/c) 2 , we get the values of da/dxp that could be compared to those presented in the 
experimental papers jlOj . In all cases, except for S + production, we overestimate the 
experimental data at small xp. 

In the region xp > 0.3-0.4, practically all the experimental xf distributions measured 
by the WA89 Collaboration [TO] are wider than the corresponding QGSM predictions ob- 
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Figure 9: The xp distributions of secondary A (left panel) and S~ (right panel) pro- 
duced in and in H~Cu interactions at 345 GeV/c, together with the corresponding 
standard QGSM calculations (solid curves) and with the QGSM calculations obtained 
by including an additional polynomial factor into the ds fragmentation functions of 
eqs. (A.20) and (A.21). 



tained by using the standard diquark fragmentation functions given by Reggeon counting 
rules. To solve this steady disagreement, we have included in the diquark fragmentation 
functions (mainly in ds diquark fragmentation) one additional polynomial factor which 
takes the spectra up at xp > 0.3-0.4. Another problem faced when trying to describe 
the experimental data of ref. [10J is connected with the behaviour at small xp. The 
peaks present in the theoretical estimations and shown in Fig. |8] are the natural result 
of dividing the rather flat function xe • ;fr by xe- The experimental behaviours of da/dx 
correspond, on the contrary, to deep minima in • ^. 

Nevertheless, the most disturbing question we face when comparing the QGSM pre- 
dictions for the Xe spectra of secondary hyperons comes from the fact that the exper- 
imental Xf distributions of secondary S~ are very different to the corresponding xp 
spectra of secondary protons in pA collisions. Since the standard QGSM predictions 
have always provided a good agreement with the experimental data on the xp spectra 
of secondary non-strange secondaries, from the theoretical point of view it would be 
puzzling if the disagreement of the corresponding spectra of strange secondaries with 
the experimental data would be confirmed. On top of that, should the theoretical pre- 
dictions for the S~ case be modified to agree with the currently available experimental 
data, this would mean significant changes in the standard QGSM diquark distributions 
and/or in the (is-diquark fragmentation function to £~. However, one actually has not 
any apparent theoretical support for those changes, and that could even result in the 
violation of the SU(3) flavour symmetry. Let's then see whether the experimental data 
on the Xp spectra for secondary £~ will be confirmed or corrected in the future. 
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Appendix 

The diquark and quark distribution functions in S~-baryon for a diagram with n cut 
Pomerons have been parametrized as follows: 

u dd (x,n) = C dd -x a *~ 2aB .(l- x y a - +m \ (A.l) 

u ds ( x , n) = C ds ■ x <*r- 2 <*b+(<*r(o)-<* v (o)) . ( X _ ? (A.2) 

u d (x, n) = C d ■ x~ aR ■ (1 - x )««- 2 «s+(««(o)-^(o))+m 1 ? ( A 3 ) 

u s (x, n) = C s - x~ a * ■ (1 - x ) aR - 2aB+m2 , (A.4) 

U sea (x, n) = C a - X~ aR ■ (1 - x )*K-^B+n-l+k _ 

The values of mi and m-i in Eqs. (A.1)-(A.4) are determined from momentum conserva- 
tion: 

(x) d + (x)ds + 2 • (n - 1) • (x) sea = 1, (A.6) 
(x) s + (x) dd + 2 ■ (n - 1) • (x) sea = 1 , (A.7) 



where 



and 



(x)i = x ■ uAx) ■ dx , (A. 8) 

Jo 



f 

Jo 



Ui(x, n) ■ dx — 1, i = q,qq,sea. (A. 9) 

/o 

For the values = 0.5, = 0, «b = —0.5, and k = 0, we obtain mi = m 2 = |(n — 1). 

The fragmentation functions of quarks and diquarks used for the description of 
strange baryon inclusive spectra are presented here. 

For quarks one has: 





(z) 


= ■ 


(1- 


_ z jX+a R -2-a B . ( 1 + fll ^m) ? 


(A.10) 




(z) 


= a A ■ 


;i - 


- z 'j x + a n- 2a B+ 2 ( 1 ^ a R)+' 2Aa . n _|_ aiZ ain ) 


(A.ll) 




(z) 


— G d 


(z) 


= ag± • (1 - z )^H-2a B +Aa . ^ + ^am) ? 


(A.12) 


GT 


(z) 


= Gf 


(z) 


= ag± • (1 - 2 ) A +^- 2a s+2(l-a fl )+Aa . ^ + ^am) ? 


(A.13) 


GT 


(z) 


= Gf 


(z) 


= ^ • GUz) , 


(A.14) 






a A 


gT 


(*) 


= aj- 


;i- 


_ z )A+a fl -2a fl +Aa.( 1+fli/1 „) ? 


(A.15) 


gT 


(z) 


= a= • 


;i - 


- z 'j x + a n- 2a B+^ Aa + 2 ( 1 - a R) . _|_ aiZ a ^^ 


(A.16) 


G° 


(z) 


= a= ■ 


;i - 




(A.17) 
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where Aa = o>r — a^, a\ = 12, a ln = 2, and A = 2a'- < pi >= 0.5. The fragmentation 
functions which are not presented here can be found in ref. [2D]. The values of the 
parameters a&, ag, ag, and a n can be obtained from a^ through quark combinatorics 
(see below). 

The QGSM diquark fragmentation functions corresponding to antibaryons production 
through the diagram shown in Fig. 2 are the following: 

G~L = G^ = G^ = a K -(l-z)W"«-^ +2 < 1 - a °»** t (A.18) 
GZu — Gud — Gdd — a s • (1 - z ) Aa ■ G^ u , (A. 19) 

Gi = a A ■ (1 - z )^R-2a B )+2(l-a B )+2Aa ? (A.20) 

G% = ^-(l-^-Gi, (A.21) 

a A 

For the diquark fragmentation functions to baryon production, and as it was men- 
tioned in the main text, they have more complicated forms than the quark fragmentation 
functions, and they contain two different contributions. The first one corresponds to the 
central production of one BB pair, and it is accounted for in Eq. (4) by fragmentation 
functions with the form: 

GL = Gt d = G K dd = a A ■ (1 - z) X-a R+i (l-a B)+ Aa ^ ^ 22) 

r Y,+ _ r T,~ _ r Z- _ r Y>+ _ r Y>+ _ r Y>- _ /-, \A-« fl +4(l-a B )+A« 

_ ^dd — ^uu — ^dd — ^ud — ^ud ~ "£± ' \ L z ) > 

(A.23) 

= ?*.(l-z)*"-G£ l , (A.24) 
a A 

a n ■ (1 - z ) a -^+ 4 < 1 -^)+ 3Aq , (A.25) 

-a fl +4(l-a fl )+2Aa 2g) 

(A.27) 



G a 


— G^d 




Gdd 


" uu 


— ^ud 




^dd ~ 


U ds 


= «A • 


(1 




U ds 


— U ds 




a A 


Gds 


a E 

a A 


(1 


-*) 


U ds 


Oq _ 

as 


(1 


-*) 



yA 

T ds 5 



Aa " Gj s , (A.28) 

A-a J j+4(l-Q S )+4Aa 

In these expressions we have used different parameterizations for the diquark frag- 
mentation functions to antibaryon production and to central baryon production, since 
antibaryons are produced one cut kink higher in the multiperipheral chain than the 
companion baryon (see Figure [2]). The corresponding expressions for proton and A pro- 
duction in pp collision were first given in ref. [2J . 

The second contribution in the diquark fragmentation functions to baryon production 
comes from the direct fragmentation of the initial baryon into the secondary one with 
conservation of SJ, shown in Figs. [3j These contributions are determined by the following 
fragmentation functions: 

G K dd = a N -z^.\v^e-(l-zf + v^z^.(l-z) 
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:i-z) 
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Act 



X-z 



,2+Act 
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'X-z 
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+ ^-^ +Aa -(l-.))+^-^-(l-^ 
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E+ 
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ud 



— ^dd 

+ Vm 

- ^ud 

+ «r 



■ z 
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{l-zf + v 



,2-/3 



2.5-/3 



^ds 



^ds 



Gdd 



G 



s- 



: CLj\i • Z 

2 ~ X-z 

p 







• e 
■fl 



(l-z) 2 .(l-z) 



A a 



P 



vf-eXl 



{l-z) 



2+Aa 



+ vi 



= ajv • z 
+ vf- 

= ^ud 



2-/3 



z) 1+Aa + v ds 



\2+Aa 



Z 



2.5-/3 



'0 



2-/3+Aa 



a+Ac 



1-2 

2-/3 _ f-| _ 



+ V ds ■ Z 



2.5-/3 



0>n ■ Z 



P 



G 



ds 



a/v • z 



P 



v, 







+ V. 



G: 



G n 



n- 

uu 



ds 
°ud 



z 



2.5-/3 



• e 
(1 



■ e ■ (1 - 2) 2 + -z 

.2-/3 



^ - a - *)] ■ a - *: 



1+Aa 



^dd 



\2+3Aa 



(7s 



QiV • z? • K • e ■ (1 - z) 1+2Aa + v ds ■ z 2 ^] -(l-z) 



1+2Aq 



(A.30) 

(A.31) 

(A.32) 

(A.33) 
(A.34) 

(A.35) 

(A.36) 
2Aa (,A.37) 

(A.38) 

(A.39) 
(A.40) 



The third term in eqs. (1A.30|) and (1A.31|) describes the contribution of the leading 
resonance E*~ and its subsequent decay into A + tt to the dd and ds fragmentation 
functions into A. The factors (1 — z) in Eq. flA.301) and (1 — |) in Eq. (IA.31j) account for 
the suppression of this contribution. 

The different probabilities for the SJ without valence quarks, Vq , the SJ with one 
valence quark, v^, and the SJ with two valence quarks, v^ q , to go into the secondary 



baryon B were deduced through quark combinatorics [261 EZ] • Here we assume that the 
strange quark suppression is common to the three diagrams shown in Fig. [3j and thus, 
e.g. for the fragmentation of the SJ without valence quarks into different baryons one 
gets: 



p : n : A + £ : S° : E~ : fi" = 4L 3 : 4L 3 : 12L 2 S : 3LS 2 : 3LS 2 : S 3 



(A.41) 



where the ratio S/L determines the strange suppression factor, and 2L + S = 1. In the 
numerical calculations we have used S/L = 0.32. 

Following the standard treatment in which the S° are included into A, and to dis- 
criminate A from the charged S, we use the empirical rule: S + + X - = 0.6 • A [27J. 

The values of Vq, v^, and used in Eq. (9) are presented in Table 1. 
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n 


A + S° 


£+ 




^0 




Q~ 


^0 


4L 3 


4L 3 


7.5L 2 5 


(9/4)L 2 5 


(9/4)L 2 5 


3L5 2 


3L5 2 


5 3 


Vu 


3L 2 


L 2 


(5/2)LS 


(3/2)L5 




5 2 






V<1 


L 2 


3L 2 


(5/2)L5 




(3/2)L5 




5 2 




V s 


_ 


_ 


(5/2)L 2 


(3/4)L 2 


(3/4)L 2 


2L5 


2L5 


5 3 


^uu 


2L 




(1/4)5 


(3/4)5 










Vud 


L 


L 


5 












Vdd 




2L 


(1/4)5 




(3/4)5 








Vus 






(5/4)L 


(3/4)L 




5 






Vds 






(5/4)L 




(3/4)L 




5 




Vss 












L 


L 


5 



Table 1: The values of parameters t> 4 in Eq. (9) obtained from quark combinatorics. 



Finally, we obtain from Eq. (IA.41j) the relations among the values of the corresponding 
parameters a# in the fragmentation functions needed in Eq. (4) for the production of 
different BB pairs: 



1 : y/(15/B) ■ (S/L) : ^(9/16) • (S/L) : ^(3/4) ■ (S/L) 2 
: v/(3/4) • (S/L) 2 : ^(1/4) ■ (S/Lf . (A.42) 
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